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ABSTRACT 

We explore the degree of magnetization at the jet base of M87 by using the observational data of 
the event horizon telescope (EHT) at 230 GHz obtained by Doeleman et al. By utilizing the method 
in Kino et al., we derive the energy densities of magnetic fields {Ub) and electrons and positrons {U±) 
in the compact region detected by EHT (EHT-region) with its full-width-half-maximum size 40 /ras. 
First, we assume that an optically-thick region for synchrotron self absorption (SSA) exists in the 
EHT-region. Then, we find that the SSA-thick region should not be too large not to overproduce the 
Poynting power at the EHT-region. The allowed ranges of the angular size and the magnetic field 
strength of the SSA-thick region are 21 //as < 0thick <26.3 //as and 50 G < i?tot < 124 G, respectively. 
Correspondingly Ub ^ U± is realized in this case. We further examine the composition of plasma and 
energy density of protons by utilizing the Faraday rotation measurement (RM) at 230 GHz obtained 
by Kuo et al. Then, we find that Ub ^ U± + Up still holds in the SSA-thick region. Second, we 
examine the case when EHT-region is fully SSA-thin. Then we find that Ub ^ U± still holds unless 
protons are relativistic. Thus, we conclude that magnetically driven jet scenario in M87 is viable 
in terms of energetics close to ISCO scale unless the EHT-region is fully SSA-thin and relativistic 
protons dominated. 

Subject headings: galaxies: active — galaxies: jets — radio continuum: galaxies —black hole physics 
—radiation mechanisms: non-thermal 


1. INTRODUCTION 

Elucidating the formation mechanism of relativistic 
jets in active galactic nuclei (AGNs) is one of the long¬ 
standing challenges in astrophysics. Although magneti¬ 
cally driven jet and wind models are widely discussed in 
the literatures (e.g., Okamoto 1974; Blandford & Znajek 
1977; Blandford & Payne 1982; Ghiueh et al. 1991; Li et 
al. 1992; Uchida 1997; Okamoto 1999; Koide et al. 2002; 
Tomimatsu & Takahashi 2003; Vlahakis & Konigl 2003; 
McKinney and Gammie 2004; Krolik et al. 2005; McK¬ 
inney 2006; Komissarov et al. 2007; Komissarov et al. 
2009; Tchekhovskoy et al. 2011; Toma & Takahara 2013; 
Nakamura and Asada 2013; McKinney et al. 2013), the 
actual value of the strength of magnetic field (B) at the 
base of the jet is still an open problem. In order to test 
magnetic jet paradigm, it is most essential to clarify the 
energy density of magnetic fields {Ub = and 

that of particles at the upstream end of the jet where 
Btot is the strength of total magnetic fields. 

Recently, short-millimeter radio observations at 
1.3 mm (equivalent to the frequency 230 GHz) have been 
performed against the nearby giant radio galaxy M87. 
M87 is located at a distance of Dl = 16.7 Mpc (Jordan 
et al. 2005; Blakeslee et al. 2009), hosts one of the most 
massive super massive black hole M, = (3 — 6) x 10® Mq 
(e.g., Macchetto et al. 1997; Gebhardt and Thomas 2009; 
Walsh et al. 2013) and thus M87 is known as the best tar¬ 
get for studying the upstream end of the jet (e.g., Junor 
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et al. 1999; Hada et al. 2011). The Schwarzschild radius 
is Rs = 2GM,/c® Ri 2 X 10^® cm for the central black 
hole with M, = 6 x 10® Mq where G is the gravitational 
constant and c is the speed of light. This corresponds to 
the angular size of ^ 7 //as. Hereafter, we set this mass 
as the fiducial one. The Event Horizon Telescope (EHT) 
composed of stations in Hawaii and the western United 
States has detected a compact region at the base of the 
M87 jet at 230 GHz with its size 40 //as (Doeleman et al. 
2012). Furthermore, Kuo et al. (2014) obtained the first 
constraint on the Faraday rotation measure {RM) for 
M87 using the submillimeter array (SMA) at 230 GHz. 

Short milli-meter VLBI observations of EHT at 
230 GHz (equivalent to 1.3 mm) is crucially beneficial in 
order to minimize the blending effect of sub-structures 
below the spatial resolutions of telescopes. Histori¬ 
cally, single dish observations of AGN jets at centi-meter 
waveband (with arc-minute spatial resolution) revealed 
that their spectra are flat at cm waveband (Owen et 
al. 1978). Marscher (1977) suggested the importance 
of VLBI observations for distinguishing various possi¬ 
ble explanations for the observed flatness. Gotten et 
al. (1980) conducted VLBI observations at cm waveband 
and found that the flat spectrum results from a blend¬ 
ing effect of sub-structures with milli-arcsecond (mas) 
scale. This was a significant forward step. However, 
subsequent VLBI observations have revealed that such 
mas scale components still have sub-structures when ob¬ 
serving them at higher spatial resolution (i.e., shorter 
wavelength). This is a vicious-circle between telescopes’ 
spatial-resolutions and sizes of sub-structures. In the 
case of M87, we finally start to overcome this problem 
since the spatial-resolution of EHT almost reaches one 
of the fundamental scales, i.e., ISGO (Innermost Stable 


2 


Kino et al. 


Circular Orbit) scale (Doeleman et al. 2012). Hence, in 
this work, we will assume the ISCO radius (i?isco) as 
the minimum size of the jet nozzle. 

Motivated by the significant observational progresses 
by EHT, we explore the magnetization degree {U±/Ub) 
in the core of M87 seen at 230 GHz. We note that Doele¬ 
man et al. (2012) did not derive i?tot and U±/Ub, and we 
will estimate them at the EHT-region for the first time. 
In the theoretical point of view, we have developed the 
methodology for the estimation of U±/Ub and Btot in 
Kino et al. (2014) (hereafter K14), and it is also appli¬ 
cable to 230 GHz. In K14, we estimated U±/Ub and 
Btot in the radio core at 43 GHz with 0 fwhm = 110 /ras 
and 0.7 Jy. We obtained the tight constraint of field 
strength (1 G < Btot < 15 G), but the resultant en¬ 
ergetics are consistent with either the 17±-dominated or 
Ds-dominated (1 x 10“® < U±/Ub < 6 x 10^). The ra¬ 
dio core at 230 GHz with 0fwhm = 40 fias directly cor¬ 
responds to the upstream end of the M87 jet and then 
it would give tightest constraints for testing the mag¬ 
netic jet paradigm. The goal of this work is applying the 
method of K14 to the EHT-detected region and exploring 
its properties. 

Opacity of EHT-region against SSA is critically impor¬ 
tant. We should emphasize that it is not clear that EHT- 
region is SSA-thick or not because short-mm VLBI ob¬ 
servations are conducted only at 230 GHz and therefore 
it is not possible to obtain spectral informations at the 
moment. Intriguingly, Rioja & Dodson (2011) detect the 
core shift between 43 and 86 GHz (in Eigure 5 in their pa¬ 
per) , which means that the radio core at 86 GHz contains 
the SSA-thick region. With the aid of interferometry 
observations, we can also infer the turnover frequency. 
The fluxes measured by IRAM at 89 GHz (Despringre et 
al. 1996) and SMA at 230 GHz (Tan et al. 2008) also 
seem to indicate that the radio core is SSA-thick above 
89 GHz although the data are not obtained simultane¬ 
ously (see also Abdo et al. 2009). The sub-mm spectrum 
obtained by ALMA also shows the spectral break above 
^ 100 GHz (Doi et al. 2013). Therefore, we may infer 
that SSA turnover frequency for the EHT-region is above 
^ 100 GHz. As a working hypothesis, we firstly assume 
that the EHT-region includes the SSA-thick region and 
apply the method of K14. We will also discuss the fully 
SSA-thin case in §6. 

The layout of this paper is as follows. In §2, we briefly 
review the method of K14. In §3, we apply the method to 
the EHT-region. In §4, the resultant U±/Ub and Btot are 
presented. In §5, we further discuss constraints on the 
proton component. In §6, we discuss the fully SSA-thin 
case. In §7, we summarize the results and give account 
of important future work to be pursued. In this work, 
we define the radio spectral index a as (x 

2. METHOD 

Eollowing K14, here we briefly review the method for 
constraining magnetic field and relativistic electrons in 
radio cores. 

2.1. Basic assumptions 

Eirst of all, we show main assumptions in this work. 

• We assume that the emission region is spherical 
with its radius R which is defined as 2R = OoSsBa 


where 0obs, Da = DlI{ 1 + z)^ and Dl are the 
observed angular diameter of the emission region, 
the angular diameter distance and the luminosity 
distance, respectively. This is justified by the fol¬ 
lowing observational suggestion. In the EHT ob¬ 
servation of M87 in 2012, Akiyama et al. (2015) 
measures the closure phase of M87 among the three 
stations (SMA, GARMA, and SMT). The closure 
phase is the sum of visibility phases on a triangle of 
three stations (e.g., Thompson et al. 2001; Lu et al. 
2012). Akiyama et al. (2015) shows that the mea¬ 
sured closure phases are close to zero (< ±20°) for 
the structure detected in Doeleman et al. (2012), 
which is naturally explained by a symmetric emis¬ 
sion region and disfavors significantly asymmetric 
one. 

• We do not include the GR effect for simplicity. The 
full GR ray-tracing and radiative transfer may be 
essential for reproducing detailed shape of black 
hole shadows (e.g., Ealcke et al. 2000; Takahashi et 
al 2004; Broderick & Loeb 2009; Nagakura & Taka¬ 
hashi 2010; Dexter et al. 2012; Lu et al. 2014). 
However, current EHT can only detect flux from 
a bright region via visibility amplitude and spatial 
structure can be constrained only by closure-phases 
(e.g., Doeleman et al. 2009). Although the pre¬ 
dicted black hole shadow images in details seem 
diverse, the size of the bright region is roughly 
comparable to ISGO scale (e.g.. Fish et al. 2013 
for review). Therefore, we do not include the GR 
effect but explore a fairly wide allowed range for 
the bright region size 0thick, he., from ~ i?isco to 
^ 2i?isco (see sec. 5). 

2.2. General consideration 

Given the SSA turnover frequency (z^ssa) and the angu¬ 
lar diameter size of the emission region at zzssa, one can 
uniquely determine Btot and K± where K± is the nor¬ 
malization factor of relativistic (non-thermal) electrons 
and positrons (e.g., Kellerman & Poliny-Toth 1969; Bur- 
bidge et al. 1974; Jones et al. 1974a, 1974b; Blandford 
& Rees 1978; Marscher 1987). Recently, K14 points out 
that the observing frequency is identical to I'ssa when 
we can identify the SSA-thick surface at observing fre¬ 
quency. 

As a first step, we assume that the EHT-region is a one- 
zone sphere with isotropic magnetic field (Utot) and par¬ 
ticle distributions in the present work. Locally, we denote 
(7 ?_l, local = Btot sin a) as the magnetic field strength per¬ 
pendicular to the direction of electron motion (Ginzburg 
& Syrovatskii 1965, hereafter GS65) where a is the pitch 
angle between the vectors of electron velocity and the 
magnetic field (e.g., Rybicki & Lightman 1979). Then, 
we can obtain pitch-angle averaged local defined as 
B± as follows: 

DL = ( 1 ) 

since B\ = (l/47r) f BtotSin^adfl = 2Btotl3. (This is 
a slightly different definition of B± in K14. The corre¬ 
sponding slight changes of numerical factors are summa¬ 
rized in Appendix.) Since we assume isotropic field, here¬ 
after we choose B± direction to the line of sight (LOS). 
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The number density distribution of relativistic elec¬ 
trons and positrons n±(e±) is defined as (e.g., Eq.3.26 
in GS65) 

n±{e±)dt± = K±t^de± (e±,min < e± < e± ,max); (2) 

where e± = j±meC^, p = 2a -|- 1, e±,min = 7±,minTOeC^, 
and e±,max = Ti.maxWeC^ are the electron energy, the 
spectral index, minimum energy, and maximum energy 
of relativistic (non-thermal) electrons and positrons, re¬ 
spectively. Although electrons and positrons may have 
different heating/accerelation process in e~ jjp mixed 
plasma (e.g., Hoshino and Arons 1991), here we assume 
that minimum energies of electrons and positrons are 
same for simplicity. By evaluating the emission at the 
synchrotron self-absorption frequency, we obtain 


2. The minimum Lorentz factor of relativistic elec¬ 
trons and positrons (yj-min) should be smaller than 
the ones radiating the observed synchrotron emis¬ 
sion (t'syn.obs), for example 230 GHz. Otherwise, 
we would not be able to observe synchrotron emis¬ 
sion at the corresponding frequency. This is gener¬ 
ally given by 

'^syiijObs ^ 1-2 X 10 m i , ’ C^) 

i + z 

These relations significantly constrain on the al¬ 
lowed values of 7±,min and Htot- 

In the next section, we will add another constraint con¬ 
dition (i.e., minimums size limit). 


B, = 




obs \ '■ 

cm) 


e. 


obs 


1 mas 


5, 


^ssa,obs 

1 Jy 


1 + z 


(3) 


where b(p) is tabulated in Marscher (1983), Hirotani 
(2005), and K14. The term K± is given by 


K±=k{p) 


Da 
1 Gpc 


-1 


X 


5, 


^ssa,obs 

1 Jy 



/ ^ssa,obs ^ / ^obs \ 

V 1 GHz / \ 1 mas / 

(5 

1 + z) 


-2p-5 


(4) 


where k[p) is tabulated in K14. The cgs units of K± and 
k[p) depend on p: ergP“^cm“^. It is useful to show the 
explicit expression of the ratio U±/Ub as follows: 

U± Wn / Da VVj^ssa.obs\-2P-13 

Ub 362(p) (p - 2) V1 Gpc/ V 1 GHz ) 

\ 1 mas / \lJy/ \l-|-zy 

(for p > 2). (5) 


From this, we find that zzssa.obs and dobs have the same 
dependence on p. Using this relation, we can estimate 
U±/Ub without minimum energy (equipartition B field) 
assumption. It is clear that the measurement of 9oh& is 
crucial for determining U±/Ub- 
We further impose two general constraint conditions. 


1. Time-averaged total power of the jet (Tjet) esti¬ 
mated by jet dynamics at large-scale should not be 
exceeded by the one at the jet base 


Tjet max [TpQy, , 

4vr 9 9 

i± = yrWct/i, 

JiTT 

ipoy = yr^/gA^ct/s, (6) 


3. APPLICATION TO THE EHT-REGION 
Here we apply the method to the EHT-region in M87. 

3.1. On basic physical quantities 

Here we list the basic physical quantities of the M87 
jet. 


• The total jet power Ljet can be estimated by con¬ 
sidering jet dynamics at well-studied bright knots 
(such as knots A, D and HST-1) located at kpc 
scale (e.g., Bicknell & Begelman 1996; Owen et al. 
2000; Stawarz et al. 2006). Based on the literatures 
on these studies, here we adopt 

1 X 10^"^ erg s-^ < Ljet < 5 x 10^^ erg s-\ (8) 

(see also Rieger & Aharonian 2012 for review). We 
note that Young et al. (2002) indicates Ljet 
3 X 10^^ erg s“^ based on the X-ray bubble struc¬ 
ture which is significantly smaller than the afore¬ 
mentioned estimate. The smallness of Ljet esti¬ 
mated by Young et al. (2002) could be attributed 
to a combination of intermittency of the jet and an 
averaging of Ljet on a long time scale of X-ray cav¬ 
ity age. In this work, we do not utilize this small 
Ljet- 

• We would assume that the bulk speed of the jet is 
in non-relativistic regime at the jet at the EHT re¬ 
gion since both theory and observations currently 
tend to indicate slow and gradual acceleration so 
that the flow reaches the relativistic speed around 
103“^ Rs (McKinney 2006; Asada & Nakamura 
2014; Hada et al. 2014). The brightness tempera¬ 
ture of the 230 GHz radio core is below the critical 
temperature ^ 10 ^^ k limited by inverse-Gompton 
catastrophe process (Kellermann & Pauliny-Toth 
1969). When the 230 GHz emission originates 
from the SSA-thick plasma, the characteristic elec¬ 
tron temperature is comparable to Tf, (e.g., Loeb & 
Waxman 2007) and Tf, at 230 GHz is in relativistic 
regime. Therefore, we set 


where L±, Lpoy, T, and fdc are, respectively, elec¬ 
tron/positron kinetic power, Poynting power, bulk 
Lorentz factor, and bulk speed of the jet at the 
EHT region. Note that Ub, U±, R are directly 
constrained by VLBI observations. 


E/? — Csound^ 




(9) 


where Csound is the sound speed of relativistic mat¬ 
ter. This will be used in Eq. m as r^/3 = 1/V3. 
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• Last, we summarize three differences between this 
work and Doeleman et al. (2012) in terms of the 
assumptions on basic physical quantities. In this 
work, we attempt to reduce assumptions and treat 
the EHT-region in a more general way. (1) In 
Doeleman et al. (2012), they assume that the EHT- 
region size is identical to the ISCO size itself which 
reflects the degree of the black hole spin. In this 
work, we do not use this assumption. (2) Doeleman 
et al. (2012) seems to focus on the SSA-thin case. 
In this work, we will investigate both SSA-thick 
and SSA-thin cases. (3) Doeleman et al. (2012) 
seems to assume ^fwhm as the physical size of the 
EHT-region. In this work, we take into account a 
deviation factor between 0 fwhm and its physical 
size (e.g., Marscher 1983). 


3 . 2 . Dijficulties for SSA-thick one-zone model 

First, we estimate the magnetic field strength in the 
EHT region by assuming that all of the EHT-region with 
^FWHM = 40 ^as is fully SSA-thick. The field strength 
of EHT-region is estimated as 


^ssa,obs \ ^ 


Htot = 3.4 X 10^ G i — I 
V230 GHz/ 


e. 


obs 


72 /ras 


^t/ssa,Qbs 

1-0 Jy 


-2 


1 + 2 


( 10 ) 


Marscher (1983) pointed out VLBI measured 0 fwhm is 
connected with true angular size 0obs by the relation 
^obs ~ 1-86 *fwhm for partially resolved sources, (see also 
Krichbaum et al. 2006; Loeb and Waxman 2007). Tak¬ 
ing such deviation into account, we examine the case of 
72 /ras = 1.8 x 40 jias for the estimate of B-field strength. 

What happens with this field strength? 


3 . 2 . 1 . Too large Poynting Power 

A severe problem arises if Btot ~ 3 x 10^ G is realized. 
Since we assume nearly isotropic random field which can 
be supported by the low linear polarization degree at 
230 Ghz (Kuo et al. 2014), the corresponding Poynting 
power is given by 

/^poy = l-5 X lO"^’^ erg s“^ 

Btot y / 2R 
300 Gy 1^1-8x1016 cm 

Here we adopt 2R = 1.8 x lO^® cm = 1.8 x 40 fias x 
16.7 Mpc. When total power of the jet (i.e., sum of 
kinetic and Poynting ones) is conserved along the jet at a 
large-scale, then this is too large compared with the jet’s 
mean kinetic power inferred from its large-scale dynamics 
a fewxlO^'i erg s“i (e.g., Rieger & Aharonian 2012 for 
review) We emphasize that a constraint on Btot by Lpoy 
is almost model-independent. 

If we allow some kind of fast magnetic reconnection 
processes (e.g., Kirk & Skj^raasen 2003; Bessho & Bhat- 
tacharjee 2007; Takamoto et al. 2012; Bessho & Bhat- 
tacharjee 2012) in order to dissipate magnetic fields at 
the EHT region, then fast and large variabilities would 
be naturally expected. However, there is no observa¬ 
tional support for such variabilities. Therefore, it seems 
difficult to realize too large Btot at EHT region. 




3 . 2 . 2 . Too fast synchrotron cooling 

Once we obtain a typical value of Btot, then we can es¬ 
timate a typical synchrotron cooling timescale. It signif¬ 
icantly characterizes observational behavior of the EHT 
region. The synchrotron cooling timescale is correspond¬ 
ingly 

1 X 10-day y^)-’. (12) 

This is much shorter than day-scale although the flux 
at 230 GHz measured by the EHT remains constant 
during subsequent three days (see the Supplementary 
Material of Doeleman et al. 2012). Then, a difficulty 
arises due to this short t±,syn- The 230 GHz radio emit¬ 
ting electrons are in the so-called fast cooling regime 
(Sari et al. 1998) in which injected electrons instan¬ 
taneously cool down by synchrotron cooling. Hence, 
a slight change/fluctuation of B field strength instan¬ 
taneously (on timescale t±,syn) is reflected on the syn¬ 
chrotron flux at EHT region. Hence, for realizing the 
observed constant flux, a constant plasma supply of Btot 
and K± with very small fluctuation is required to avoid 
rapid variability/decrease of synchrotron flux. On the 
other hand, when the magnetic fields are not that large, 
t±,syn can become longer than day scale. Then, we can 
avoid rapid variability/decrease of synchrotron flux with¬ 
out imposing very small fluctuation of Btot and K± in 
the bulk flow. Since some fine-tuning of Btot and K± in¬ 
jection may be able to adjust the observed constant flux 
density at the EHT region, the too-fast-cooling problem 
may be less severe than the aforementioned problem on 
too-large-Lpoy. But it is natural to suppose that smaller 
Btot realizes in the EHT region to avoid fine tuning of 
injection quantities. 

3 . 3 . Two-zone model 
3 . 3 . 1 . Basic idea 

The difficulty of too-large-Lpoy can be resolved if the 
EHT-regions is composed of SSA-thick and SSA-thin re¬ 
gions and the angular-size of SSA-thick region (0thick) is 
more compact than Oohs, i-C., 

^obs ^ ^thick- (i^8) 

We show an illustration of our scenario in Figure [TJ In 
this solution, most of the correlated flux density detected 
by EHT is attributed to the emission from the SSA-thin 
region. Since Vssa. of the SSA-thin region is by dehnition 
smaller than 230 GHz, the magnetic field must be signif¬ 
icantly smaller because the held strength is proportional 
to Because of this reason, we regard the SSA-thick 
region as the main carrier of the Poynting power. 

Here, we assume the ISCO radius for non-rotating 
black hole (i?isco = 6GM,/c^ = 3i?s = Oi^coDl) as the 
minimum size of SSA-thick region. This corresponds to 
the angular size 21 ^as. Indeed, theoretical works (Brod¬ 
erick & Loeb 2009; Lu et al. 2014) comparing EHT ob¬ 
servations and jet models also indicate the model images 
with short-mm bright region ^ with their size comparable 
to ISCO. Therefore, we examine the range of SSA-thick 
region 6»thick > 21 /ras. 

^ Conventionally, such regions are sometimes called as hot spots 
in literatures (e.g., Lu et al. 2014 and reference therein) 
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3.3.2. Gaussian fitting with two-components 

In Figure m we estimate the correlated flux density of 
this SSA-thick region based on the EHT data. The ob¬ 
served flux data plotted as a function of baseline length 
are adopted from Doeleman et al. (2012). The black 
solid curve is the best-fit circular gaussian model by 
Doeleman et al. (2012). The red solid curve is the 
best-ht model. The red dashed and dot-dashed curves 
represent the SSA-thick and the SSA-thin components, 
respectively. 

Below, we explain the details of the gaussian fitting. 
To determine the correlated flux density for the com¬ 
pact SSA-thick region with its lower limit size 0fwhm = 
21 yLtas/1.8 = 11.1 /ras, we conduct the two-component 
(SSA-thick and thin components) gaussian fitting to the 
EHT-data. First, we obtain the upper limit of the 
correlated flux density for the SSA-thick component as 
Si, = 0.27 Jy. Next, we perform the the two-component 
(SSA-thick and thin components) gaussian fitting by fix¬ 
ing 0FWHM = 21 ^as/1.8 = 11.1 /ras and S^, = 0.27 Jy. 
Then, we obtain the corresponding size and flux of 
the extended SSA-thin component S^, = 0.75 Jy and 
^FWHM = 60 /ias. 

4. RESULTS 


region can be estimated as 


n = 


27rkiy^„„ „K„(0thick/2)2 


syn,obs 

!2 X 10^° K 


obs 


^thick 

0.27 Jyy ^^21 ^as 


-2 


(15) 


where JZsyn.obs = 230 GHz. This value is comparable with 
the Tb at 86 GHz estimated by Lee (2013). 

Last, it is worth to add one thing. According to the 
equation of state in relativistic temperature regime (e.g., 
Ghandrasekhar 1967; Kato, Fukue & Mineshige 1998), 
we obtain ® 


TijSsa 


o kTb 

nieC^ 


10 


^thick 

21 ^as 


-2 


(16) 


where we use the fact that Ti^minWeC^ can be identical 
to the average energy of electrons and positrons since p is 
steeper than 2. The obtained 7±,ssa tends to be smaller 
than the minimum Lorentz factor obtained by Eq. m 
by a factor of a few. While we may use y± ssa as 7±,min, 
we conservatively use the condition Eq. (Ildl) taking some 
uncertainty of numerical factor in Eq. m into account. 


Here, we limit on Btot, ^thick, U±/Ub, in the EHT- 
region without assuming plasma composition. The crit¬ 
ical value 7±.min is derived by the combination of the 
jet power limit (Eq. |6]) and synchrotron emission limit 
(Eq. JT])). By eliminating Btot, we obtain 

7±,min<1.2 X 10^ 

‘2R ( Ljet 

1.8 X 10^® cm J \ 5 X 10^^ erg s“^ J ’ 

(14) 

where t'ssa = 230 GHz is used. Since 7±,min has R de¬ 
pendence, larger R allows slightly larger 7±,min- 

In Eigure O we show the value of log(C/±/t/_B) in 
the allowed ranges of 7±,min and Btot with Ljet = 5 x 
10^"^ erg s“^ and p = 3.0. It is essential to note that the 
maximum value of Btot is determined by the condition 
Bpoy < Bjet while the minimum value of Btot is gov¬ 
erned by the condition of 0thick > Bisco/B*l ~ 21 pa,s. 
The right side of the allowed region is determined by the 
I'syn.obs limit shown in Eq. ©. Note that the maximum 
value of 0thick = 26.3 pas is smaller than 40 pas. This 
suggests that the EHT-region has a more compact SSA- 
thick component in it. Interestingly, overall SSA-thick 
region satisfies Ub ^ U±. If protons do not contribute 
to jet energetics, then this result supports the magneti¬ 
cally driven jet scenario. In Table[T]we show the resultant 
allowed values. Summing up, we find that (1) the allowed 
^thick satisfies 21 pas < 0thick < 26.3 pas, and that (2) 
the allowed fields strength is 50 G < Btot < 124 G. 

When we choose a smaller Bj, the upper limit of ^thick 
and Btot becomes smaller according to Eqs. jS]) and ([6]). 
When Bj = 1 X 10"^'^ erg s“^, the allowed Btot and ^thick 
are 50 G < Btot < 65 G and 21 pas < 0thick < 22.4 pas. 
In this case, the allowed regions of Btot and dthick are 
very narrow. 

We add to note a short comment on brightness tem¬ 
perature. The brightness temperature of the SSA-thick 


5. CONSTRAINTS ON PROTON COMPONENT 

In § [5l we investigate constraint on the energy density 
of protons (Bp) by using Faraday RM measured by Kuo 
et al. (2014). From the measured RM, we will constrain 
the number density of protons (rip). Then, we exam¬ 
ine Up. The degree of proton contribution in energetics 
has a significant influence over relativistic jet formation 
(e.g., Begelman, Blandford & Rees 1984; Reynolds et al. 
1996). 

5.1. Further Assumptions 

To discuss the proton contribution, we need to add 
some further assumptions. Although the observed ra¬ 
dio emissions warrant the existence of relativistic e~ 
population, it is not clear about the origin of relativistic 
e~ je'^ which radiate radio emissions at 230 GHz. There 
are several possibilities for its origin. Relativistic pro¬ 
tons may play an important role for heating/acceleration 
of positrons via resonance process with relativistic pro¬ 
tons in shocked regions (e.g., Hoshino and Arons 1991), 
while direct pair injection (Iwamoto & Takahara 2002; 
Asano & Takahara 2009), and/or relativistic neutron in¬ 
jection (Toma & Takahara 2012) processes may also work 
at the jet formation regions. It is beyond the scope of 
this work to clarify the origin of relativistic e~ popu¬ 
lation and their relation with proton component. In this 
section, we simply assume the existence of protons and 
generally define the average energy of these protons as 

Cp. 

As mentioned in the Introduction, Kuo et al. (2014) 
obtained the first constraint on RM for M87 using SMA 
at 230 GHz. Although it is not clear how much fraction of 
linearly polarized emission comes from the EHT-region, 
it is worth to extend the method used in the previous 

® When magnetic fields are uniform, the numerical factor at 
the right-hand side in Eq. m is smaller than this case because 
of fewer degree of freedom for electrons/positrons (see Slysh 1992; 
Tsang and Kirk 2007). 
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sections by including RM constraint and apply to the 
present case of 230 GHz core of M87. The degree of LP 
^ 1 % at 230 GHz detected by Kuo et al (2014) is sig- 
nihcantly smaller than the value when fully-ordered field 
(i.e., typically ^ 70 % for SSA-thin case and ^ 16 % 
for SSA-thick case, see Pacholczyk 1970). Hence, the 
assumption of isotropic H-fields in this work looks rea¬ 
sonable to some extent. On the other hand, only ordered 
magnetic fields aligned to the line of sight (Hlos) con¬ 
tribute to the RM. Hereafter, we conservatively assume 
Btot > 73los- 


In sub-section 5.4, we will constrain proton contribu¬ 
tions in the case of i?tot ~ Hlos in Eq ([19]). At the 
moment, this is the only case which we can deal with 
within this simple framework. 

5.3. Plasma composition and /p-coupling rate 

To further constrain physical properties at the jet base, 
here we introduce the basic plasma properties and define 
general notations. The number densities of protons (rip) 
positrons (n+), and electrons (n_) are, respectively, de¬ 
fined as follows: 


5.2. RM limit 

Here we introduce a new constraint using RM obser¬ 
vation data. This RM is important for estimating the 
kinetic power of protons (Tp) because RM can constrain 
the proton number density. Generally speaking, an ob¬ 
served rotation measure (AMobs) consists of two parts, 
i.e., RM by internal (jet) (AMjet) and RM by external 
(foreground) matter (RMext)- Therefore, the i?Mobs can 
be decomposed into 

RAfobs = RMjet + RMext- (17) 


Up = rjn ^, 

n+ = (l —i?)n_ (0 < 7 ? < 1 ), 

np = He-p ~ 9 _ „ „ _ 1 '^±7±finim (^l) 

Z IIP 1 

where 17 is a free parameter describing the proton-loading 
in the jet. Here we use the charge neutrality condition 
in the jet. It is convenient to define further quantities: 

n± = n- + = (2 — r])n-, 


Basically, it is difficult to decouple RM^^t and RMext and 
obtain RMjet. However, it may be possible to discuss 
an upper limit of |i?A/jet| with some reasonable assump¬ 
tions. When the observed RM (RMohs) is comparable to 
RMext, then we obtain 


RMobs ~ RMext, iRMjetl < RMohs- (18) 

Indeed, foreground Faraday screen in close vicinity of 
jets seems to well explain observed RMohs for radio-loud 
AGNs (e.g., Zavala & Taylor 2004) . The explicit form 
of RMjet the rotation measure for relativistic plasma is 
given as 




27rm2c^ 


dlBhos^p 


log 7 ±,ti 


±,min 


< 5.36 X lO^Htot^p 


log 7±,i 
27:“ 


±,min 


R 


10 ^® cm 


rad 


(19) 


where we set f dl ^ 2R since the region is assumed as 
uniform. From this, we see that Faraday rotation is 
strongly suppressed in relativistic plasma (Jones & Odell 
1977; Quataert & Gruzinov 2000; Broderick & McKinney 
2010). Note that RM only include ionic plasma contri¬ 
bution, and does not include the electron/positron pair 
plasma. It is because electron and positron have the same 
mass but have opposite (i.e., minus and plus) charges and 
then the net Faraday rotation by them is cancelled out. 
Qualitatively saying, the mixture of pair plasma (i.e., 
77 < 1 ) effectively reduce the value of RMjet- 
Regarding i?M-limit of M87, Kuo et al. (2014) has 
measured \RMobs\ ~ (3.4— 7.5) x 10 ® rad m“^ and they 
assume RMohs ~ RMgxt- Following Kuo et al. (2014), 
we also assume RMohs ~ RMext- Then, the RM-limit 
can be written as 


I RMjet I < I X 10® rad m-^. (20) 

Note that the above constraint only gives the upper 
limit of Tip. Therefore, the finite value of |RMjet| does 
not exclude the plasma composition of pure plasma. 


where n± and n^-p are, the number density of electrons 
and positrons, and that of proton-associated electrons, 
respectively. The case of 77 = 0 corresponds to pure 
plasma while 77 = 1 corresponds to the pure e~ jp 
plasma. Next, it is important to clarify energy balances 
between electrons and protons. It is useful to introduce 
the parameter defining the average energy ratio between 
protons and electrons as C as 

TT7 

e± = Cep, (— < C < 1), (23) 

TRp 

where e± is the average energy of relativistic e^. The 
case C = 1 can be realized for equipartition between 
electrons, positrons and protons via effective /p cou¬ 
pling while C = TTie/mp means inefficient e^ jp coupling 
for example through randomization of bulk kinetic en¬ 
ergy of the jet flow (e.g.. Kino et al. 2012 and reference 
therein). Since we focus on the case of p > 2 suggested in 
M87 (Doi et al. 2013), relativistic electrons at minimum 
Lorentz factors characterize the total energetics. Here, 
e± Ri 7±,minTOeC^ Can be estimated as 0.5 MeV < e± < 
50 MeV together with 1 < 7 ±,min < 100 based on the ob¬ 
tained 7 ±,inin- Then, the case C = 1 corresponds to that 
of non-relativistic protons (0.5 MeV < < 50 MeV) 

while the case C = Tne/mp coincides with that of rela¬ 
tivistic protons (1 GeV < Cp < 100 GeV). 

In general, L is decomposed to 

Tjet — Lp Tpoy, 

L± = L_+L+, (24) 

where L±, L_, L+, Lp, and Tpoy are, the powers of the 
sum of electrons and positrons, electrons, positrons, pro¬ 
tons, and magnetic fields respectively. For convenience, 
we define for Lp = L± and it is given by 

2 C 

Lp > L± for 
Lp = L± for 
Lp < L± for 


V > rieq, 
rj = 77eq,and 

77 < 77eq. (25) 
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Since we set 


U±Kie±n±, (26) 

Up^CpThp — p, 

Lp/L± = Up/U± = ??/[(2 — 77 )C] holds. Finally, time- 
averaged total power of the jet (ijet) can be generalized 
as follows: 


Ljet > max 


L 


poy, 






(28) 


Given the two model parameters rj and C, we obtain Up. 


5.4. Limits on Btot, ^'thick, U±/Ub, and Up/Us 
Here, we give limits on Btot, ^'thick, U±/Ub, and Up/UB 
in the EHT-region for e~ je^jp mixed plasma. As for 
plasma properties, the following four cases with proton 
loaded plasma can be considered, i.e., relativistic protons 
with e“/p-dominated composition, relativistic protons 
with e^-dominated composition, non-relativistic protons 
with e“/p-dominated composition, and non-relativistic 
protons with e^-dominated composition. 


5.4.1. The case for relativistic protons (C, = 


Here we consider the case for relativistic protons (^ = 
nie/mp). In Figure |4l we show a typical example of 
”e“/p-dominated” case with 77 = 0.99. In this case, 
we obtain r^gq = 1.09 x 10“^. Since we consider ”e“/p- 
dominated” composition, the upper limit of RM signifi¬ 
cantly constrains smaller 7 ±,min according to Eq. (IT^ . In 
this case, Ub U± still holds as smaller 7 ± min region is 
excluded by the RM constraint. In Table [21 we summa¬ 
rize the resultant allowed physical quantities in this case. 
The maximum value of Btot is determined by the condi¬ 
tion Lpoy < Ljet while the minimum value of Btot is gov¬ 
erned by the condition that 0 thick > 7 ?isco/.Dl ~ 21 pas. 
In the limit of inefficient /p coupling, minimum energy 
of electrons/positrons are smaller than that of protons 
by a factor of mg/mp (i.e., e±,min = («^e/mp)ep,min)- 
Therefore, L± decreases and Lp tends to dominate over 
L±. The energetics constraint in this case is given by 


Ljet > max 


L 


poy 5 


/1 I_ n 

' 2—n nie 



In the case of e^-dominated composition with smaller 
rj also leads to the same Btot and U±/Ub. In the same 
way as shown above, the maximum and minimum val¬ 
ues of Btot are determined by the jet power limit and 
minimum size limit at the EHT-region. However, Up is 
much smaller than Ub simply because of the paucity of 
the relativistic proton component. 


5.4.2. The cases for non-relativistic protons {( = IJ 


Next, let us consider the case of non-relativistic protons 
(( = 1). When non-relativistic protons are loaded, the 
corresponding energetic condition can be given by Ljet > 


max 


^’( 1 + 2 ^) 


L± 


Since the protons are non- 


relativistic, the effect of proton loading is quite small in 
terms of energetics. The coefficient resides in a narrow 
range 1 < (1 -I- 77/(2 — 77 )) < 3/2. Note that RM strongly 
depends on 77 while RM is independent of C- 
The ”e“/p-dominated” case results in similar values 
of Btot and U±/Ub to those shown in Table [21 because 


the maximum and minimum values of Btot are also de¬ 
termined by the jet power limit and minimum size limit 
at the EHT-region. The contribution of protons are only 
Up = U±/2. So, it does not give any significant effects 
on energetics. 

Einally, we comment on the ”e^-dominated” case. The 
main difference between the ”e“/p-dominated” and ”e^- 
dominated” cases is n^-p- Since the number density of 
e^-pairs does not contribute to RM, the constraint of 
RM becomes weaker when n^-p becomes smaller. It 
leads to wider allowed region for smaller 7 ± and smaller 
B region. Therefore, the maximum value of allowed 
U±/Ub for the ” e^-dominated” case becomes larger than 
that for ”e“/p-dominated” case. However, this only 
changes the allowed 7 ± within a factor of ~ 10 and it 
does not give a large impact on energetics. 


6. FULLY SSA-THIN CASE 

It is worthwhile to examine a case of fully SSA-thin 
model for EHT-region since the indication of Vssa, > 
100 GHz by interferometry observations does not nec¬ 
essarily mean that Vss& is larger than 230 GHz. We can 
safely regard the SSA frequency as 43 GHz < r'ssa < 
230 GHz where the lower limit is warranted by the de¬ 
tection of core-shift at 43 GHz in Hada et al. (2011). 

In Figure O we show a schematic draw of the syn¬ 
chrotron spectrum when the EHT-region is SSA-thin at 
230 GHz (solid line). The upper limit of the flux den¬ 
sity at 43 GHz of the 230 GHz core is estimated as 
0.09 Jy = 0.7 Jy x (40/110)^ based on the VLBA mea¬ 
surements of the radio core flux and size by Hada et al. 
(2013). The gray-colored scale shows the typical flux 
density obtained by SMA and GARMA. Interferometric 
observation shows some variability at 230 GHz (Akiyama 
et al. 2015). We define this as Fupper and we assume that 
Fupper is the upper limit of the flux density in overall fre¬ 
quency range of 43 GHz < r'ssa < 230 GHz. First, from 
the EHT data, we can estimate a possible lower limit of 

^ssa 

7...a>230 GHz X f 

V 5j./l Jy 

-160 GHz (for a = 2.5). (29) 



Note that, the choice of a = 3.0 leads to zzssa 170 GHz. 
Second, from the VLBA data, we can estimate a possible 
upper limit of T/gga as 


^ssa < 43 GHz X 
-160 GHz. 


/ f upper/2. 3 Jy 

V S'^/0.09Jy 


2/5 


(30) 


Allowing some flux measurement errors, somehow we can 
have consistent case around i^ssa ^ 160 GHz with a — 

2.5. 

Then, let us discuss on physical quantities in this case. 
From Eq. O U±/Ub oc Vssa^~^^■ Therefore, in this 
case, the ratio would be typically larger by a factor of 
(160/230)“^® — 6.9 X 10^ (for p = 2.5) than that for 
the SSA-thick case. However, this does not change the 
result of U± <C Ub since C/± <C Gb < 10“"^ in any 
cases with the SSA-thick core existing. Hence, we can 
conclude that even for fully SSA-thin EHT-region case. 
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U± Ub holds in order not to overproduce fluxes be¬ 
tween 43 GHz < i/gsa < 230 GHz. 

However, a critical difference appears for the compar¬ 
ison between Up and Ub- In the case of relativistic pro¬ 
tons with ”e“/p-dominated composition, Up > Ub can 
be realized for a certain range of t'ssa- From the Table [2] 
we know the values of Up/Us when i/ssa = 230 GHz. By 
multiplying the factor of ^ 200 — 400 , and the maximum 
value reaches Up/Us > 1 at Vssa ^ 160 GHz. 

7. SUMMARY 

We have explored the magnetization degree of the jet 
base of M87 based on the observational data of the EHT 
obtained by Doeleman et al. (2012). Following the 
method in K14, we estimate the energy densities of mag¬ 
netic fields {Ub) and electrons and positrons {U±) in the 
region detected by EHT (EHT-region) with its FWHM 
size 40 /ras. Imposing basic energetics of the M87 jet, the 
constraints from EHT observational data, and the min¬ 
imum size of the SSA-thick region as the ISCO radius, 
we find the followings. 

• First, we adopt the assumption that the EHT- 
region contains an SSA-thick region. Then, the 
co-existence of SSA-thick and SSA-thin regions is 
required for the EHT-region not to overproduce 
Lpoy. The angular size of the SSA-thick region 
is limited as 21 /ras < ^thick < 25.5 ^as, while 
that of the SSA-thin region should be 40 ^as to 
explain the EHT data. The derived flux density of 
the SSA-thick region is about 0.27 Jy. The allowed 
magnetic-fields strength in the SSA-thick region is 
58 G < Htot < 127 G. In terms of energetics, 
Ub ^ U± is realized at the overall SSA-thick re¬ 
gion. If protons do not dominantly contribute to jet 
energetics, then this result supports the magnetic- 
driven jet scenario at the SSA-thick region. 

We further examine the following four cases 
for electron/positron/proton {e~ je^jp) mixed 
plasma; non-relativistic protons with e~ jp- 
dominated composition, non-relativistic protons 
with e^-dominated composition, relativistic pro¬ 
tons with e“/p-dominated composition, and rel¬ 
ativistic protons with e^-dominated composition, 
together with the assumption that RM detected 
by SMA (Kuo et al. 2014) gives an upper limit 
of RM of the EHT-region. Although RM limit 
can give tighter constraints on allowed 7 ±, it does 
not change the results significantly. We find that 
Ub ^ U± always holds in any case. 

• Second, the case of completely SSA-thin {I'ssa < 
230 GHz) EHT-region is also discussed. Although 
lower Vssa can increase the ratio U±/Ub by a fac¬ 
tor of 200 — 400 than that for the SSA-thick case, 
this does not change the result of U± <C Ub since 
U±/Ub < 10“^. However, we also find that, in the 
case of relativistic protons with ”e“/p-dominated” 
composition. Up > Ub can be realized around 
I'ssa ~ 160 GHz. 

Future work and key questions are enumerated below. 

• An important future work is to confirm the exis¬ 
tence of the SSA-thick region in the EHT-region. 


If we confirm it, then we can exclude the case of 
Up/Ub >1- In the context of conhrming the ex¬ 
istence of SSA-thick region, we also add to note 
the effectiveness of inclusions of longer baselines 
even for a single frequency VLBI observation. In 
Fig [21 it is clear that the visibility amplitude of 
the SSA-thin component is much smaller than that 
of SSA-thick component above ~ 3GA at 1.3mm 
wavelength. Therefore, inclusions of baselines with 
> 3GA would be effective to distinguish the SSA- 
thick component. For example, phased ALMA plus 
SMT with an effective bandwidth of 4 GHz would 
be effective at ~ 5GA (Fig 6 in Fish et al. 2013). 
In Figure |2l we show the corresponding baseline- 
length range (the blue-shaded region). 

• Equally important future work is to observe the 
EHT-region with the spatial resolution of ^ 1 i?s 
of M87. Currently, the EHT array with 20-30 pas 
resolution at 230 and 345 GHz (e.g., Lu et al. 2014) 
is not able to reach ~ 1 i?s of M87. Ground- 
based short-mm VLBI observations are very sen¬ 
sitive to weather conditions (e.g., Thompson et 
al. 2001). To confirm our assumption that the 
minimum OthickDL is comparable to ~ i?isco or 
even smaller, space VLBI observations would be 
required in future. In the past missions and exist¬ 
ing plan of space VLBI, it was not able to reach 
the event horizon scale of M87 (e.g., Dodson et al. 
2006; Asada et al. 2009; Dodson et al. 2013) since 
target wavelength were not short enough. Thus, 
atmospheric-free space (sub-)mm VLBI observa¬ 
tion would be indispensable to reach ^ 1 of 
M87. The phased ALMA (e.g., Alef et al. 2013, 
Fish et al. 2013) will play a definitive role for 
such observations for obtaining visibilities between 
space and ground telescopes baselines. 

Honma et al. (2014) have recently proposed a new 
technique of VLBI data-analysis to obtain super¬ 
resolution images with radio interferometry using 
sparse modeling. The usage of the sparse model¬ 
ing enables us to obtain super-resolution images in 
which structure finer than the standard beam size 
can be recognized. A test simulation for imaging 
of the jet base of M87 is actually demonstrated in 
Honma et al. (2014) and the technique works well. 
Therefore, this super-resolution technique will be¬ 
come another important tool for obtaining better 
resolution images. 

• The observational result of Doeleman et al. (2012) 
does not show flux variability at 230 GHz. How¬ 
ever, total epoch-number of EHT observations is 
too scarce to confirm the absence of flux variability 
at 230 GHz all of the time. M87 might be in qui¬ 
escent state during the EHT observations in April 
2010 by chance. We also emphasize that the de¬ 
rived held strength is still > 58 G and t±,syn still 
tends to be smaller than day scale. It is also in¬ 
triguing that the same correlated hux densities in 
2009 reported by Doeleman et al. (2012) are ob¬ 
served during another EHT observation performed 
in April 2012 (Akiyama et al. 2015). This result 
is quite different from the day scale variability de- 
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tected in Sagittarius A* by the EHT observations 
(Fish et al. 2011). To search for a possible flux 
variability of M87 in more details, continuous mon¬ 
itoring by EHT would be essential. 

• Based on GRMHD model, well-ordered poloidal 
fields are dominant within the Alfven point and 
toroidal fields become dominant outside of Alfven 
point while turbulence may not grow-up yet at the 
jet base (e.g.. Spruit 2010 for review). In general, 
turbulent eddies which most probably generate tur¬ 
bulent fields are not expected before sufficient in¬ 
teractions with surrounding ambient matter (e.g., 
Mizuta et al. 2010 and reference therein). There¬ 
fore, higher LP degree is likely to be expected. 
Conservatively saying, the reason of low LP degree 
by Kuo et al. (2014) is most probably because of 
depolarization within SMA beam. At the moment, 
we are not able to rule out a possible constitution of 
RIAF emission which may also lead to low LP de¬ 
gree. If so, then studies of fundamental process for 
particle accelerations in RIAF (e.g., Hoshino 2013) 
and the effects particle escape from RIAF (Le & 
Becker 2004; Toma & Takahara 2012; Kimura et 
al. 2014) would become more important. 

• In terms of the brightness temperature of the 

230 GHz radio core of M87 Tf, ~ 2 x 

lOio k(^) seems slightly higher 

than the prediction of hot electron temperature of 
~ 10® K in RIAF flows (e.g., Manmoto et al. 1997). 
Hence, the jet emission seems to be preferred to ex¬ 
plain EHT-emission in M87 (Dexter et al. (2012), 
see Ulvestad & Ho (2001) for similar arguments). 
However, it is not conclusive because geometry near 
ISCO regions is highly uncertain in observational 
point of view. The scrutiny of the origin of the 
230 GHz emission is still a noteworthy big issue to 
explore. 

• Furtherpolarimetric observation would be required 
to examine RM properties in more details. Al¬ 
though we adopt RM values of Kuo at al. (2014), 
it is found that the observed electric vector po¬ 
sition angle (EVPA) trend does not show a suf¬ 
ficiently tight fit to A®-law. This behavior may not 
be due to the consequence of blending of multiple 
cross-polarized sub-structures with different RM 
values, but simply rather due to the non-uniformity 
between the upper and lower side bands of the 


SMA. A polarimetric observation with ALMA is 
clearly one of the promising first step to improve 
this point. Obviously, in the final stage, short-mm 
(and sub-mm) VLBI polarimetric observations are 
inevitable to avoid the contamination from the ex¬ 
tended region. 

• Degree of e^/p coupling is a critical factor for 
the results of the proton power. Theoretically, 
Hoshino and Arons (1991) found the energy trans¬ 
fer process from protons to positrons via absorp¬ 
tions of high harmonic ion cyclotron waves emitted 
by the protons. Amato and Arons (2006) indeed 
performed one-dimensional particle-in-cell (PIG) 
simulations for e“/e“''/p-mixed plasma. However, 
there are several simplifications in PIC simulations 
such as smaller irip/me ratio etc. More intensive 
investigations are awaited to clarify the degree of 
e^/pcoupling at the base of the M87 jet. 

• We make a brief comment on effects of magnetic 
field topology and anisotropy of e“/e+ in terms 
of energy distribution. If e“/e+ energy distri¬ 
bution in the EHT region is isotropic, then the 
synchrotron absorption coefficient investigated by 
GS65 is applicable and differences of field-geometry 
would not have an impact on field strength estima¬ 
tion. For example, the difference of i?tot between 
the cases of isotropic held (see Eq. (1)) and ordered 
held {B± = Btot) which is directed towards LOS is 
only a factor of a/3/2. 

However, if the e“/e^ energy distribution is highly 
anisotropic, then the well known synchrotron emis- 
sivity and self-absorption coefficient are not appli¬ 
cable. Effects of the e“/e+ anisotropy on syn¬ 
chrotron radiation are not well studied and it is 
beyond the scope of this paper. Although we 
do not have any observational suggestions of such 
anisotropy of e~ je^ energy distribution, it may be 
a new theoretical topic to be explored if observa¬ 
tional suggestions are found in the future. 
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APPENDIX 

MODIFICATION OF NUMERICAL FACTORS 

In order to obtain better accuracy calculation plus some modihcations of the dehnition of B± and relevant corrections, 
modihed numerical co-efficient of b{p) and k(p) are presented although the corrections are small. 

In K14, magnetic fields strength perpendicular to the local electron motions were not averaged over the pitch angle 
(In Eqation (1) in K14). In this work, in Eqation ([T]), we conduct the pitch-angle averaging for dehning the averaged 
magnetic fields strength perpendicular to the local electron motions. 

Synchrotron self-absorption coefficient measured in the comoving frame is given by (Eqs. 4.18 and 4.19 in GS65; Eq. 
6.53 in Rybicki & Lightman 1979) 


CXi/ 


•/Se^ / 3e 
SniTie \27rmfc^ 


Cl (p) K± B^p+ 2) i/- (p+4) /2 ^ 


(Al) 
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where the numerical coefficient ci (p) is expressed by using the gamma-functions as follows; ci(p) = r[(3p-|-2)/12]r[(3p-|- 

22)/12]. For convenience, we define . 

Optically thin synchrotron emissivity per unit frequency from uniform emitting region is given by (Eqs. 4.59 and 
4.60 in BG70; see also Eqs. 3.28, 3.31 and 3.32 in GS65) 


= 47r- 




/ 3e \ 


(p-l)/2 


(A2) 


where the numerical coefficient is C2(p) = r[(3p -I- 19)/12]r[(3p — l)/12)]r[(p -|- 5)/4)]/r[(p -|- 7)/4)]/(p -|- 1). In K14, 
Btot was wrongly written as B±. So, here we revise it and it leads to larger b{p) by the factor of vTTs. For convenience, 
we define ei, = 'iTTX 2 C 2 {p)B[yy^'^^‘^. The modified coefficient b{p) is expressed as 


b{p) = 


/47rc2A2 X 1.5^/^ 


6ciXi 


X 2 


-4 


(A3) 


In K14, the index of square bracket at the right hand side of b{k) should not be 2 but —2(typo). The expression of 
k{p) oc does not change, but the value k{p) is changed. Although the modifications of b and k in Table 1 

of K14 are straightforward based on the above explanations, we put the table [3] for convenience. 
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optical ly*thin 
region (40 pas) 


Fig. 1.— Illustration of the jet base of M87 down to the EHT-region scale. The right panel shows the actual image of M87 with VLBA at 
43 GHz adopted from Hada et al. (2013). The yellow-green circle shows the one-zone region with its diameter 110 /i.as which is investigated 
in K14. The EHT-region detected by Doeleman et al. (2012) is shown as the blue circle. Since Hada et al. (2011) indicate that the central 
engine of M87 is located at ~ 41 fi&s eastward of the radio core at 43 GHz, we put the the EHT-region around there. The left panel shows 
the illustration of internal structure inside the EHT-region. The red-colored region represents an SSA-thick compact region inside the 
SSA-thin region. The black-colored region conceptually shows a possible BH shadow image. According to the smallness of closure phase 
reported in Akiyama et al. (2015), a certain level of symmetry is kept in this picture. 



Fig. 2.— Gaussian fittings to the correlated flux density of the M87 core obtained by EHT at 230 GHz. The flux density data plotted 
as a function of baseline length are adopted from Doeleman et al. (2012). The black solid curve is the best-fit circular Gaussian model 
with Si/ = 0.98 Jy and ^fwhm = 40 /ias obtained by Doeleman et al. (2012). The red solid curve is the best-fit two-component model. 
The red dashed and dot-dashed curves represent the SSA-thick and the SSA-thin components, respectively. The SSA-thick component is 
expressed as the gaussian with ^fwhm = 21 /i.as/1.8 = 11.1 /ias and Su = 0.27 Jy. The size and the flux density of the extended SSA-thin 
component are ^fwhm = 60 /ias and Su = 0.75 Jy. The blue-shaded region represents the baseline-length range corresponding to the one 
between the Hawaii/Arizona/Galifornia and Chile. 


TABLE 1 

Results when the EHT region contains SSA-thick region 


L 

allowed Btot 

allowed 0thick 

allowed U±IUb 

[erg s 1] 

[G] 

[/ias] 


5 X 10^^ 

50 < Stot < 124 

21 < ethick < 26.3 

7.9 X 10"’' < -^ < 2.3 X 10-3 
— Ub — 
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Fig. 3.— The allowed region of 7±,min and Btot (the red cross points enclosed by the black trapezoid). The colored contour lines show the 
allowed \o^{U±/U-q). The tags \o^{U±/U-q) =-4, -4.4, -5, and -5.4 are marked as reference values. The physical quantities and parameters 
adopted are Ljet = 5 x ergs~^, and p = 3.0. The minimum is limited by t'syiijobs at 230 GHz. 



Fig. 4.— The allowed region of 7±,min and Btot when RM limit is taken into account. The physical quantities and parameters adopted 
are Ljet = 5 x 10'^'^ ergs~^, p = 3.0, p = 0.99 and = melrrip which corresponds to e“/p-dominated composition with relativistic protons. 
The tags log(C/j-/C/B) = -4.4, -4.8, -5.2, -5.4, and -5.6 are marked as reference values. 



log [frequency (GHz)] 

Fig. 5.— Schematic view of synchrotron spectrum when the EHT-region is fully SSA-thin at 230 GHz with its size and flux density 40 pas 
and 1.0 Jy (solid line). The upper limit on the flux density at 43 GHz is estimated as 0.09 Jy = 0.7 Jy x (40/110)^ based on the VLBA 
measurements of core size and flux at 43 GHz by Hada et al. (2013). The gray-colored range shows the typical flux density at 230 GHz 
obtained by SMA and CARMA (e.g., Doeleman et al. 2012; Kuo et al. 2014; Akiyama et al. 2015). 
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TABLE 2 

Results for the case of e“/p-DOMiNATED composition with relativistic protons 
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fraction 

[%] 

Li 

[erg s 1] 

allowed Btot 

[G] 

allowed 0thick 
[/ras] 

allowed U±IUb 

allowed UplUs 

0.9 

10 

5 X lO'*^ 

50 < Rtot < 124 

21 < ethick < 26.3 

7.9 X 10“^ < ^ < 1.1 X lO-"' 

— Ub — 

1.2 X 10-3 < < 0.17 

— Or — 

0.99 

1 

5 X 10^4 

50 < Rtot < 124 

21 < ethick < 26.3 

7.9 X 10"^ < ^ < 1.1 X 10-4 
— Ub — 

1.4 X 10-3 < ^ < 0.20 
— Ub — 

1 

0 

5 X 10'*'“ 

50 < Rtot < 124 

21 < ethick < 26.3 

7.9 X 10“^ < ^ < 1.1 X 10-4 
— Ub — 

1.4 X 10-3 < ^ < 0.21 
— Ub — 


TABLE 3 

Relevant coefficients for Rt and K± 


p 

Kp) 

b(p) in K14 

h{p) in Hirotani (2005) 

b{p) in Marscher (1983) 

k(p) 

2.5 

4.1 X 10-5 

3.3 X 10-5 

2.36 X 10-5 

3.6 X 10-5 

9.3 X 10-3 

3.0 

2.4 X 10-5 

1.9 X 10-5 

2.08 X 10-5 

3.8 X 10-5 

1.4 X 10-3 

3.5 

1.5 X 10-5 

1.2 X 10-5 

1.78 X 10-5 

- 

2.1 X lO-"' 








